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Magnetization Step in Spatially Distorted Heisenberg Kagome Antiferromagnets 
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Motivated by recent experiment for volborthite, a typical spin- 1/2 antiferromagnet with 
kagome lattice structure, we study magnetization process of classical Heisenberg model on 
a spatially distorted kagome lattice using Monte Carlo (MC) method. We find a distortion- 
induced magnetization step at low temperatures and low magnetic fields. Magnitude of this 
step is given by Am2 = |1 — a| /3a at zero temperature, where a denotes spatial anisotropy in 
exchange constants. The magnetization step signals a first-order transition at low temperatures, 
between two phases distinguished by distinct and well developed short-range spin correlations, 
one characterized by coplanar spin alignment with a \/3 x period, and the other by a non- 
coplanar and partially spin flopped structure. We point out the relevance of our results to the 
unconventional steps observed in volborthite. 
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In antiferromagnets on the kagome lattice, geometrical 
frustration effects are believed to suppress the conven- 
tional magnetic long-range order and induce large de- 
generacy of the ground state. It opens a possibility of 
realizing spin-liquid states, spin analogues of liquids or 
ices.^ Volborthite Cu3V207(OH)2-2H20 is known as a 
nearly ideal spin-1/2 kagome antiferromagnet.^ Measure- 
ments of the magnetization and the specific heat for this 
compound indicate absence of the conventional magnetic 
long-range order down to 50 mK.^'^ Therefore, it is pro- 
posed that this material offers an evidence of the quan- 
tum spin-liquid state in nature.^ 

Recent studies on volborthite have revealed the exis- 
tence of unconventional three steps in the magnetization 
curve. ^ These steps are not anticipated in the magneti- 
zation process of isotropic kagome Heisenberg model in 
the literature^' ^ and the origin of these steps is not un- 
derstood. A MC study of the classical Heisenberg model 
has clarified the existence of the magnetization plateau at 
one-third of the saturation magnetization (mgat/S) at fi- 
nite temperatures interpreted as the order by disorder ef- 
fect.^ Even for the quantum Heisenberg model {S = 1/2 
and 1), exact diagonalization study has shown similar 
magnetization plateau at zero temperature.^ However, 
the magnetization step has not been found at all in the 
isotropic kagome Heisenberg model. 

In this letter, to shed light on the origin of the un- 
conventional steps observed in volborthite, we calculate 
the magnetization curve in the antiferromagnetic classi- 
cal Heisenberg model on the spatially distorted kagome 
lattice by using MC method. In our viewpoint, struc- 
tural distortion, which inevitably exists in volborthite,"* 
is a key to explain the magnetization step. This struc- 
tural distortion induces a spatial anisotropy in exchange 
interactions. Therefore, as a simple model, we employ 
an anisotropic (J1-J2) kagome Heisenberg model to de- 
scribe the magnetization process observed in volborthite. 
We find distortion-induced first-order phase transitions 
at low temperatures and low magnetic fields. In the vicin- 
ity of this first-order phase transition, step-like behaviors 



of the magnetization curve are found and spin structure 
factors change drastically between coplanar and non- 
coplanar structures. We also determine the h-T phase 
diagram for the spatially anisotropic kagome Heisenberg 
model. 
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Fig. 1. (Color online) An L X L (L = 2a) distorted kagome lat- 
tice. Lattice sites labeled by A, B, and C constitute a unit cell 
(dashed-line triangle) . The number of spins is TV = 3L^ . Bonds 
along two directions have exchange constant Ji (thin blue lines), 
while the bond along the third direction has exchange constant 
J2 = aJi (thick red lines). For volborthite, it is plausible that 
a > 1. 



Our model Hamiltonian is defined in the form 

= Ji ^ 5, • + J2 ^ S, ■ S, - , (1) 

where Si denotes the Heisenberg classical spin, h de- 
notes an external magnetic field, and and (1,^)2 
denote nearest- neighbor sites on AB/AC chains and BC 
chains respectively. In the kagome lattice, unit cell con- 
tains three sites A, B and C, as shown in Fig. 1. The spa- 
tial anisotropy induces different antiferromagnetic cou- 
pling on each triangle: two Ji bonds and one J2 bond. 
For simplicity, we set Ji = 1 and J2 = aJi. We simplify 
the Hamiltonian as 
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{Sa ■ Sb + O^SB ■ Sc~^ Sc • Sa) 
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The summations are taken over upward and downward 
triangles. The lattice has no distortion when a = 1. 
There are two simple limits for this model. For a ^ 0, 
the lattice becomes a decorated square lattice, with ad- 
ditional sites at the midpoints of square lattice edges. 
In this limit, the ground state is ferrimagnetic. While 
for a — >■ cxj, the lattice becomes isolated antiferromag- 
netic chains and free spins, which is equivalent to the 
one-dimensional lattice. In volborthite, we infer that a 
is larger than unity from comparison of bond lengths. 
The magnitude of the anisotropy is estimated to be less 
than 20% from theoretical analysis.'' We consider both 
kinds of anisotropy, especially a = 1.05 and a = 0.95, 
but perform detailed analyses on the a = 1.05 case. 

We clarify thermodynamic properties of this model 
by using MC simulations. The Metropolis algorithm has 
been employed. We start the simulations from random 
initial states. Spins are updated in the sequential order 
on the lattice. We also perform one over-relaxation up- 
date per one MC step to accelerate efficient samphng.^ 
At every temperature and field, 5 x 10^ MC steps were 
taken for equilibration. These are followed by measure- 
ments done during more than 5 x 10^^ MC steps. Finally, 
the results have been averaged over more than 8 inde- 
pendent runs to estimate statistical errors. At low tem- 
peratures T < 0.01, we use exchange MC algorithm.^ 
Simulations were done for the system size up to L = 72 
{N = 15552) under the periodic boundary conditions. 
We set the Boltzmann constant fee = 1. 

In Fig. 2 (a), we show magnetization rUz parallel to 
applied fields at several choices of temperatures as a 
function of magnetic fields. For both a > 1 and a < 1, 
we have observed a magnetization plateau TOsat/3 at low 
temperatures, which have already been reported by Zhit- 
omirsky for isotropic (a — 1) kagome lattice. Due to 
the distortion, the magnetization plateau becomes much 
wider. This tendency is consistent with the results of 
the exact diagonalization for spin- 1/2 distorted kagome 
Heisenberg model. ^ 

As shown in Fig. 2 (b), we have observed step- like 
behaviors at low temperatures and low magnetic fields. 
To see whether these step-like behaviors are the first- 
order phase transitions in the thermodynamic limit, we 
calculate the energy distribution P{E) up to L = 72 
{N = 15552). If the transition is of the first-order, the 
energy distribution P{E) should be bimodal at the tran- 
sition temperature. As shown in Fig. 2 (c), we find a clear 
two-peak structure around the transition point. This re- 
sult shows an evidence for the first-order phase transi- 
tion at low temperatures. These first-order phase transi- 
tions appear only when the lattice has the anisotropics 
{a 7^ 1). The transition field increases monotonically 
with increasing temperatures. The first-order jump be- 
comes reduced with increasing temperatures and van- 
ishes at a finite-temperature critical point as we detail 
below. 

The origin of the first-order phase transitions is well 
understood by considering zero temperature limit. To an- 
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Fig. 2. (Color online) (a) Magnetization curves at T = 0.01 for 
a = 1.05 and a = 0.95. The system size is L = 36. The dashed- 
dotted line denotes one third of the saturation magnetization. 
For both anisotropies, magnetization plateaus exist at h ~ 2. (b) 
Magnetization curves at low magnetic fields. Step-like behaviors 
in the magnetization curve arc observed for both a = 1.05 and 
a = 0.95. The system size is L = 36. The step shows first-order 
transitions below the critical temperature Tc ~ 0.007. At T > Tc, 
the step vanishes and turns into a crossover. We have observed 
a magnetization hysteresis for T < 0.008. Statistical errors are 
smaller than the symbol size at T = 0.008. Magnetization curves 
in the decreasing field are shown at T < 0.006. Due to finite size 
effects, a first-order step at T = 0.004 appears to be rounded. 
The dashed line denotes the magnetization at T = given by 
"n^z = (h + 2\1 — a\)/{Qa). See the text for details, (c) Energy 
distribution P{E) at several magnetic fields for T = 0.006. The 
system size is L = 60. It shows two-peak structure around the 
transition point (sec the case of h = 0.07). 



alyze the ground-state properties, we rewrite the Hamil 
tonian as 
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where Sa = Sx/ot + S-b + Sc- For /i = 0, the en- 
ergy is minimized with Sa = 0. The ground state, 
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which satisfies the constraint S/x = 0, has an infi- 
nite degeneracy and the magnetization is zero. While 
for < h < min(2,4a — 2), the energy is minimized 
with Sa = {h/2a)ez, 5*^ = — sign(l — a), and hence 
Sl^ = Sc = [h + 2sign(l - q;)]/(4q;). Therefore, at zero 
temperature, we obtain the magnetization as 



/i + 2|l 



a 



„ . • (4) 

ba 

Equation (4) indicates that the magnetization jumps by 
infinitesimal magnetic field at T = 0. Although thermal 
fluctuations reduce the magnitude of this jump and shift 
the transition point to higher fields, this first-order phase 
transition still survives at sufficiently low temperatures 
as we have already shown by MC calculations. 

From the above analyses, it turns out that the ori- 
gin of the first-order phase transition is the switching of 
the ground state by a level crossing. When the spatial 
anisotropy exists, macroscopic degeneracy of the kagome 
antiferromagnet is partially lifted by the infinitesimal 
magnetic field, i.e., under the magnetic fields, the or- 
dered vector becomes q = (0, Qy), where the range of qy 
is given by — tt < Qy < tt. 

To identify the nature of the first-order transition, we 
calculate at T > the spin structure factors defined as 



, (5) 



where the index I denotes the positions of the three sites 
on the unit cell, and i, j, and Rtj denote the coordinate 
of the (i, j)-th unit cell on the triangular Bravais lattice. 
We use the dimensionless units in which the lattice pe- 
riod (length of the square in Fig. 1) is taken unity and 
hence the interspin distance equals 1/2. The first Bril- 
louin zone is a square bounded by — tt < Qx^ Qy t^- 

As shown in Fig. 3, the structure factors S^'^{q) dras- 
tically change when we cross the first-order transition 
point. At lower fields [h = 0.1, T = 0.01, see Fig. 3 (b)], 
S^y{q) has broad peaks at q = (±27r/3, ^27r/3), which 
corresponds to -v/S x \/i pattern with the 120° struc- 
ture in real space. At higher fields [h = 0.2, T = 0.01, 
see Fig. 3 (c)], the structure factor has sharp peaks at 
q = (0,±7r), signalizing the real space pattern shown in 
the lower panel of Fig. 3 (c). In the magnetic structures 
in Fig. 3 (c), spins on A sites align parallel to the mag- 
netic field indicating a partial spin flip under magnetic 
fields. Along a BC chain, the xy components of spins on 
the B sites are all parallel each other and are antiparal- 
lel to all the spins on the C sites. Along an AB chain, 
the spins on the B sites align alternatingly. Although the 
ground states are degenerate among any ordered vec- 
tor q = (O^Qy), thermal fluctuations appear to favor the 
q ~ (0, ±7r) configuration via the order by disorder ef- 
fect.iO'ii 

We obtain phase boundaries for a = 1.05 at lower 
temperatures and lower fields h < 2. From the loca- 
tion of the steps, we determine the first-order transi- 
tion line. Because this first-order transition does not in- 
volve any explicit symmetry breaking, its nature is es- 
sentially the same as that of gas-liquid type. Thus, the 
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Fig. 3. (Color online) (a) Magnetization curve around the first- 
order transition at T = 0.01. (b) and (c) Spin structure factors 
for a = 1.05, L = 36 (upper panel) and corresponding schematic 
spin configurations in real space (lower panel) around the first- 
order transition. In (b), h is lower than the transition field he ~ 
0.16 (h = 0.1, T = 0.01), where the structure factor shows a 
X pattern. In (c), h is higher than he (h = 0.2, T = 
0.01), where the structure factor shows a qr = (0, itTr) pattern. 
In spin structure factors, we illustrate the region —it < ija; < tt 
and —IT < < TT in the momentum space. In schematic spin 
configurations illustrated in the lower panels of (b) and (c), the 
unit cells of each pattern are shown in dotted-line. 



first-order phase transition line terminates at the criti- 
cal point and a crossover line exists beyond the critical 
point. This crossover line is characterized by the peak 
of the specific heat and the susceptibility. The univer- 
sality class of the critical point of the gas-liquid transi- 
tion belongs to that of two-dimensional Ising model. 
Therefore, along the crossover line, the susceptibility x 
has the singularity Xpoak oc \T ~ Td''^ or \h ~ hd''^, 
where 7 = 7/4 near the critical point. By using this re- 
lation, we roughly estimate location of the critical point 
as {K,T,) = [0.11(1), 0.007(1)]. 

Although continuous symmetry breaking is prohibited 
in two dimensions at T 7^ 0, a quasi-long-range order 
may occur for the case of the XY anisotropy. Because 
the effective symmetry of the Heisenberg model becomes 
XY type under the magnetic fields, the quasi-long-range 
order, i.e., Berezinsky-Kosterlitz-Thouless (BKT) tran- 
sition may occur, ^^"^^ although at T = the system 
becomes effectively disconnected BC chains. In order to 
examine the possibility of the quasi-long-range order, we 
have calculated the spin correlation function and esti- 
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mated the correlation lengths below the crossover line. 
We find that the spin correlations along the AB and BC 
chains decay exponentially for system size L > 48 even 
at the lowest accessible temperature. Therefore, we con- 
clude that no BKT phase exists under magnetic fields. 

We now examine the relevance of our results to the un- 
conventional steps observed in volborthite. In our model, 
wc find one distortion-induced magnetization step. This 
first-order transition is similar to the step-like behaviors 
in volborthite, though three steps are observed in the ex- 
periment. To identify the relationship between the exper- 
iment and the result of our model, we estimate the upper 
bound of the magnetization step. Since the distortion is 
20% at most, we obtain Aruz = (a — l)/(3a) < 1/18 at 
T = from eq. (4). At T 7^ 0, the jump occurs between 
two nonzero m^'s at h ^ and as a first look, it would 
not specify the magnetization at the step. Nevertheless 
the jump size Am^ is less than the jump at T = and 
after the jump with increasing field h is at most twice 
of Am^ at T = 0. By considering the magnitude of the 
three steps Wgat/S, 'TT-sat/6, and TOsat/45 in the experi- 
ment, if we assign our jump to one of the experimental 
jumps, it is likely that the msat/45 step corresponds to 
the present distortion-induced step. The sudden change 
in the spin structure factors from q = (±27r/3, =F27r/3) to 
q = (0,±7r) structures at the first-order transition may 
be detected in neutron scattering and nuclear magnetic 
resonance experiments. 

In the experiment, at the step does not sensitively 
depend on T, while our results imply that it should de- 
crease to niz = when T — > 0. This discrepancy may be 
due to quantum effects under which q = (±27r/3, =F27r/3) 
structure may be stabilized even at h ^ and T = 0. It 
is intriguing to identify the critical point for volborthite 
and see whether the first-order transition survives under 
quantum fluctuations. 
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Fig. 4. (Color online) (a) Temperature dependence of specific 
heat at several magnetic fields which are labeled from ft = 0.10 
to 0.45 starting from the left (a = 1.05, L = 60). (b) Phase 
diagram of distorted kagome lattice (a = 1.05). The first-order 
phase transition line terminates at the critical point (open circle) . 
Dashed line represents a crossover line, which is characterized by 
the peak of the specific heat. 



In summary, by using the MC method, we have clar- 
ified thermodynamic properties of the classical Heisen- 
berg antiferromagnets on a distorted kagome lattice un- 
der magnetic fields. The distortion-induced magnetiza- 
tion step appears at low temperatures and low magnetic 
fields. Estimating the upper bound of the step size, we 
conclude that this is the first-order transition induced by 
the spatial anisotropy and corresponds to the step at the 
lowest field observed for volborthite. The spin structure 
factor shows a sharp change at the transition which may 
be detected experimentally. 

Clarifying the origins of two other steps in the ex- 
periment remains for future studies. Other interactions 
in volborthite such as next-nearest neighbor interactions 
and Dzyaloshinsky-Moriya interactions may give rise to 
additional unconventional phase transitions. Besides, our 
results are obtained for the classical model while volbor- 
thite consists of quantum spin-1/2 spins. It is an interest- 
ing future issue whether the quantum fiuctuations sup- 
press the first-order transition observed in the classical 
system, leading to a crossover with a similarity to the 
experiment where temperature dependence of the tran- 
sition field is weak. 
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